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Abstract 

We propose a generalization of Heisenberg picture quantum mechanics in which a Lagrangian and 
Hamiltonian dynamics is formulated directly for dynamical systems on a manifold with non-commuting 
coordinates, which act as operators on an underlying Hilbert space. This is accomplished by defining the 
Lagrangian and Hamiltonian as the real part of a graded total trace over the underlying Hilbert space, 
permitting a consistent definition of the first variational derivative with respect to a general operator-valued 
coordinate. The Hamiltonian form of the equations is expressed in terms of a generalized bracket operation, 
which is conjectured to obey a Jacobi identity. The formalism permits the natural implementation of gauge 
invariancc under operator-valued gauge transformations. When an operator Hamiltonian exists as well as a 
total trace Hamiltonian, as is generally the case in complex quantum mechanics, one can make an operator 
gauge transformation from the Heisenberg to the Schrodinger picture. When applied to complex quantum 
mechanical systems with one bosonic or fermionic degree of freedom, the formalism gives the usual operator 
equations of motion, with the canonical commutation relations emerging as constraints associated with the 
operator gauge invariance. More generally, our methods permit the formulation of quaternionic quantum 
field theories with operator-valued gauge invariance, in which we conjecture that the operator constraints 
act as a generalization of the usual canonical commutators. 



1. Introduction 

The search for unification of the laws of physics has been closely intertwined with the discovery of wider 
classes of symmetries of the fundamental equations. Thus, the unification of electricity and magnetism in the 
Maxwell equations directly relates to both Abelian gauge invariance and relativistic kinematics; the advent of 
the standard model has as its underpinning the widening of the gauge principle from Abelian to non-Abelian 
groups. Current attempts at further unification are based largely on the use of supersymmetries relating 
bosonic to fermionic degrees of freedom. 

In this paper we explore another possible direction for broadening the underlying symmetries, by gen- 
eralizing from c-number to unitary operator gauge invariance of the fundamental equations. This approach 
grew out of an investigation [1] of quaternionic quantum mechanics and quaternionic quantum field theory, 
but the basic concepts can be understood without going through most of the details of quaternionic Hilbert 
space, and are presented here in a self-contained fashion. We begin, in Sec. 2, by developing an operator 
dynamics based on the idea of a total trace Lagrangian and Hamiltonian, and show how this can be used 
to formulate operator field equations which are covariant under operator-valued gauge transformations. We 
then apply the formalism in Sec. 3 to complex quantum mechanics, and show that it contains, and gener- 
alizes, the usual canonical formalism. In Sec. 4, we explain a few simple facts about quaternionic Hilbert 
space, and then apply the total trace formalism to construct operator gauge invariant quaternionic field 
theories, the properties of which are discussed. Some brief concluding remarks are given in Sec. 5. Further 
details of the material presented in this paper can be found in the final two chapters of [1]. 
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2. Total trace Lagrangian formulation of quantum 
dynamics and operator gauge invariance 

The primary tool for achieving invariance under operator-valued gauge transformations is the concept 
of a total trace Lagrangian, which will be developed in general form in this section. We begin by introducing 
an underlying complex or quaternionic Hilbert space, Vu, which we assume to be the direct sum 

Vh = V+(B (la) 

of a Hilbert space of bosonic states and a Hilbert space of fermionic states. Following Witten [2] 
we define an operator (—1)^ which counts fermion number modulo two, that is, (—1)^ has eigenvalue +1 
on all states in V^, and has eigenvalue —1 on all states in V^. Using this operator, we then define a trace 
operation Tr O for a general operator O as follows, 

Tr = Rerr(-l)^C> = Re^ {n\{-lfO\n) , (lb) 

n 

with Re the real part, with {|n)} any complete set of states, and with Tr the usual operator trace or diagonal 
sum. The operation Tr has the following useful properties: 

(i) If = £)- is fermionic, then Tr = 0, since 

Tr = Re Tr[{~lfO-] = Re Tr[0-{-lf] = -Re Tr[{-lfO-] = -Tr . (Ic) 

(ii) If O = 0+ is bosonic, and 0+ = C(i)0(2), then 0(i) and 0(2) are either both bosonic or both fermionic, 
and we have 

Tr 0(i)0(2) = Re rr[(-l)^0(i)0(2)] = Re rr[0(2)(-l)^0(i)] 

(Id) 

= ±Re Tr[(-l)^0(2)0(i)] = ±1^ 0(2)0(i) , 
with the + sign holding when and 0(2) are both bosonic, and the — sign holding when 0(1) and 

0(2) Stie both fermionic. 

(iii) If O = —O^ is anti-self-adjoint, then 

Tr O = Re Tr[{-lfO] = Re Tr[{-lfO]^ = Re Tr[{-lfO^ = -Tr O , (le) 
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and Tr O vanishes. Correspondingly, if O is self-adjoint, then Tr O agrees with Tr {—1)^0, which is 
already real. 

(iv) If Ti"^^ OrSqr = for arbitrary independent operator variations 6qr, then each Or must vanish, while 

r 

if OrSqr — for operator variations 5qr restricted to be of either bosonic or fermionic type, then 

r 

the part of Or of the same type must vanish. The first statement follows from 

Tr Ordqr = Re i {n\{-lfOr\m){m\6qr\n)] ; (If) 

r \n,m,r / 

choosing {m\Sqr\n) = (n|(— 1)^0^1"^) (with the bar denoting complex or quaternion conjugation, the 
latter as defined in Sec. 4 below) gives 

Tr^ a5g.= ^ |(n|(-l)^a|m)|2 , (Ig) 

r n,Tn,r 

which can vanish only if (—1)^0^ — 0, which implies Or = 0. The second statement follows by noting 
that when 5qr is of bosonic or fermionic type, then property (i) implies that 

TV ^ OrSqr = Tr ^ O^'^Sqr = Tr ^ O^'^ [5qr + Sqi"^] , (Ih) 

r r r 

with oi^^ the part of Or of the same type as 5qr and Sqr°^ an arbitrary variation of the opposite type 
as Sqr- But Sqr + Sqr°^ is an unrestricted variation, so the first statement of property (iv) then implies 

O^'^ = 0. 

Let now {qr{t)} be a set of time-dependent quantum variables, which act as operators on the underlying 
Hilbert space, with each individual qr of either bosonic or fermionic type, and let {qr{t)} be their time 
derivatives. We introduce an operator Lagrangian L which is a polynomial function (or more generally, a 
Laurent series expandable function) of the variables {qr} and {qr}, 

L = L[{qr}, {qr}] , (2a) 

and we define the total trace Lagrangian L by 

UUr}, {Qr}] = Tr L[{qr}, {qr}] , (2b) 
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and the total trace action S by 

/oc 
dtL . (2c) 
-oo 

Because of property (i) of Tr, any fermionic part of L is automatically projected to zero, so there is no loss 
in generality in assuming that L is bosonic. Similarly, by property (iii) of Tr, any anti-self-adjoint part of 
L is automatically projected to zero, so we lose no generality by further specifying that L is self-adjoint. 

Let us now examine the consequences of requiring the total trace action to be stationary under arbitrary 
operator variations of the {qr}, subject to the restriction that 5qr be of the same bosonic or fermionic type 
as Qr- When we vary a given variable Qr, the variation of L consists of a sum of terms of the form 

OLSqrOn , (3a) 

with Ol,r operators appearing respectively on the left (L) and right [R] of Sq,-, which in general do not 
commute with 5qr. Inside the operation Tr, we can cyclically permute the factors in Eq. (3a) to get 

Tr OLSqrOn = ±Tr OnOLSqr , (3b) 

with the +(— ) sign corresponding, as in property (ii), to whether On is of bosonic (fermionic) type. Re- 
ordering all terms with the general form of Eq. (3a) this way, we are able to identify a well-defined operator 
6Ij/6qr, of the same bosonic or fermionic type as g^) which obeys 

AT 

5L = Tr — 6qr . (3c) 
dqr 

Similarly, varying one of the qr, we identify a well-defined operator SJ-i/Sqr, again of the same type as qr, 
which obeys 

AT 

6l, = Tr — 6qr . (3d) 
dqr 

When qr is of bosonic type, the order of the factors within Tr in Eqs. (3c, d) is irrelevant, but when qr is 
of fermionic type the factor ordering is significant, since by property (ii) of Eq. (Id), a minus sign appears 
when the order of two factors of fermionic type is reversed. In many applications, some of the qr are 
either self-adjoint or anti-self-adjoint in character. If 6qr is further restricted to have the same adjointness 
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character as qr, then only the parts of Sh/Sqr and SL/Sqr which have the same (opposite) adjointness 
character as a bosonic (fermionic) qr are well-defined. It will be assumed henceforth that for those qr with 
definite adjointness character, the variational derivatives SL/Sqr and SL/Sqr denote the operators of the 
same (opposite) adjointness character as a bosonic (fermionic) qr which obey Eqs. (3c, d). We note, finally, 
that the procedure just described cannot be extended to higher order variational derivatives. Since Sh/Sqr 
is already an operator, a further variation will involve a sum of terms of the form of Eq. (3a), in which the 
Sqr (or Sqr) factors are sandwiched between operators on left and right with which they do not commute. 
Without the trace, there is now no way to combine the terms in the sum into a single expression with 
infinitesimals on the right, and hence there is no definition of second variational derivatives analogous to 
Eqs. (3c,d). 

Let us now impose an action principle, by requiring 

SS = (4a) 
under arbitrary same-type operator variations. Varying all the arguments qr and qr of L, we get 

/oo /■CXJ / CT srr \ 

dt6L= / dtTrV l — Sqr + ^dqr] , (4b) 
J-oo ^ \Hr 5qr J 



which by an integration by parts becomes 

r°° „ TAT, rl /AT.M 

5qr . (4c) 



^ Sqr 



+ 



r dtTr V I"—- — (—\ 

J-oo * ^ [s^ di [s^J 



So if we take the variations 6qr to vanish at ±oo, requiring (5S = gives, by property (iv), the operator 
equations of motion (the generalized Euler-Lagrange equations) 



6qr dt \dqr 



. (4d) 



Corresponding to the Lagrangian form of the equations in Eq. (4d), we can set up a Hamiltonian form by 
the usual method of making a Legendre transformation. Defining the momentum operator Pr conjugate to 

Qr by 



we define the total trace Hamiltonian H by 

H = Tr ^ Pr-Qr - L . (5b) 

r 

We then have, under general same-type operator variations, 

(5H = Tr {SprQr +PrSqr) ~ X/ ( ^^"^ A~ ^'^'^ ) ' ^^^^ 

which substituting Eqs. (5a) and (4d), and using property (ii), becomes 

SH = Tr {±qr6pr — Pr^Qr) , (5d) 

r 

with the +(— ) sign chosen according as whether Qr is of bosonic (fermionic) type. Equation (5d) shows that 
H is a total trace functional of the operators {qr} and {pr}, 

H = U[{qr}, {pr}] , (5e) 

with the operator variational derivatives 

su su 

T— = -Pr , — = ±qr . (5f) 

Oqr OPr 

As in the case of the Lagrangian variations, when qr has a definite adjointness character, the variations 
SH/Sqr and SH/Spr denote the operators obeying Eq. (5d) which have the same (opposite) adjointness 
character as a bosonic (fermionic) qr- We note, finally, that with pr defined as in Eq. (5a), if the Euler- 

Lagrange equations are satisfied but arbitrary variations dqr are permitted at t = ±(X), then Eq. (4c) implies 
that 



(5g) 



r 

This formula and the others involving S have obvious generalizations when the time interval (— oo, oo) is 
replaced by (Ti,T2), for arbitrary finite Ti^2- 

Continuing in analogy with the standard Hamiltonian formalism, let A[{gr}, {Pr}] and B[{gr}, {Pr}] be 
any two total trace functionals of the operator arguments {qr} and {pr}, and let us define their generalized 
Poisson bracket 

{A,B} = Tr V ± — ^ ^ , 6a 

' \Sqr 6pr 6qr 6prJ 



with the +(— ) sign again corresponding to bosonic (fermionic) . Then for a general total trace functional 
A[{q'r}, {Pr}] we have 

{A,H} = T.^(±)(-----j=TVi: (^a. + ^..j=^A, (6b) 

and since by construction the generalized Poisson bracket is anti-symmetric in its arguments, 

{A,B} = -{B,A} , (6c) 

it follows that the time derivative of H vanishes, 

|h = {H,H} = 0. (6d) 

An important, and perhaps difficult, open question concerning the generalized Poisson bracket of Eq. (6a) 
is whether it satisfies a Jacobi identity. That is, let A[{gr}, {Pr}]) B[{gr},{pr}] and C[{g'r}, {Pr}] be any 
three total trace functional of the operator arguments {qr\, {Pr\, and let us define the bracket 

[A, B, C] = {A, {B, C}} + {C, {A, B}} + {B, {C, A}} , (6e) 

which is totally antisymmetric in A, B and C. On the basis of a number of examples, we conjecture that 

[A,B,C] = 0, (6f) 

but we do not have a proof. It clearly is an important issue to find either a proof or a counter-example. If 
the Jacobi identity is satisfied, then using Eq. (6b), the bracket {A, B} of any two constants of the motion 
A and B is itself a constant of the motion. 

We now have all the ingredients needed to give a generalized version of Heisenberg picture quantum 
mechanics. States are described by fixed vectors \b) € and |/) e V^, and so the inner product geometry 
specified by the set of all inner products {{b\h')} and {(/|/')} is automatically time independent. The time 
dependence of the operators {q^} and {p^} is completely specified by Eq. (5f), giving these operators at all 
times once their form is specified at some initial time (say t = 0). The most general observable O will be a 
self-adjoint polynomial function (or Laurent expandable function) of {qr}, {Pr} and the time t, 

= 0[{qr},{pr}A, (7a) 



and its time dependence is determined by using the Leibnitz product rule and Eq. (5f). The expectation of 
O in any state |6) or |/) can be rewritten as a total trace functional according to 



{O) 



{{h\0\h)=TrPbO 

\ {f\0\f) = -Tr PfO , 

Pt = \h){h\, Pf = \f){f\. (7b) 



This permits us to apply Eq. (6b), as generalized to the case in which A has an explicit time dependence, 

giving 

I {0) = {dO/dt) + {{0),H}. (7c) 
Transition probabilities can also be reexpressed as total trace functional, 

\{b\b')\^ = Tr PbPy = Tr Pb'Pb , \{f\f)f = -TrPfPf, = -TrPf,Pf , (7d) 

and are time independent by virtue of the time independence of the projectors Pb,Pf,... . 

After this rather lengthy excursion into total trace Lagrangians and Hamiltonians, we are ready to 
introduce the concept of operator gauge invariance. In its simplest form, an operator gauge transformation 
consists of a transformation on the operators Qr of the form 

qr UrqrUl + ^qr[Ur\ , (8a) 
with each Uj. a unitary operator of bosonic type, 

UrU:^ = UlUr = 1 , [(-1)^, Ur]=0, (8b) 

and with Ag^fC^r] an inhomogeneous term calculable in terms of the operator Uj.- The Lagrangian L of 
Eq. (2a) cannot in general be constructed to be invariant under the transformation of Eq. (8a), but we will 
find that we can readily construct Lagrangians L in the form 

L = Y, , (8c) 

r 

which transform under Eq. (8a) as 

L^Y^ UrLrUl . (8d) 

r 
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As a consequence, although the operator Lagrangian L is not invariant, the total trace Lagrangian 



L = TrL (8e) 

is invariant under Eq. (8a), 

L ^ L , (8f) 

by virtue of the properties of Ur in Eq. (8b) together with the cyclic invariance of the trace. We will also 
employ a second form of operator gauge transformation, in which the variables Qr divide into three groups, 
the operators Qr in the first group transforming as in Eq. (8a), those in the second group transforming as 

Qr ^ U^QrK'' + Aqr[U^] , (9a) 

and those in the third group transforming as 

Qr ^ UrQrK^ , (9b) 

with Ur, U'r two independent unitary operators of bosonic type. We will now find Lagrangians L in the form 

L = ^ (L, + 4) , (9c) 

r 

which transform under Eqs. (8a) and (9a,b) as 

r 

Again, while the operator Lagrangian L of Eq. (9c) is not operator gauge invariant, the corresponding total 
trace Lagrangian L is operator gauge invariant. 

We now make a number of remarks concerning the structure and properties of operator gauge invariant 
total trace Lagrangians. 

(1) We have written the equations of this section with r a discrete index, but in many of the applications 
described in the next two sections, r will be replaced by a continuum coordinate x. 
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(2) Varying Eqs. (8a), (9a) and (9b) with respect to Qr, the inhomogeneous term drops out, and we get 
respectively 

6qr UrSqrU}. (first group) , 
Sqr —>■ U'^SqrU'^ (second group) , 

5qr Ur5qrU!} (third group) . (9e) 

Hence when L is operator gauge invariant, the Eulerian derivative appearing in Eq. (4c) transforms as 

Er — > UrErU^ (first group) , 

Er — > U^ErU^^ (second group) , 

Er UlErUl (third group) , 
(5L d 



(9f) 



6qr dt Sqr 

and the Euler-Lagrange equations Er = are operator gauge covariant. 
(3) A total trace version of the familiar Noether theorem of classical mechanics can be derived, as follows. 
Let 6A{t) be an infinitesimal operator parameterizing a set of operator variations 6qr of the variables 
qr, and let us assume that SL only involves 6A{t) and SA{t), but not SA{t) or higher time derivatives 
(the Lagrangians studied in the next two sections all have this feature). Then we have 



SL = Tr{^^SA+^-^6A 
SA 5A 



5S = TV ^ M 
5A 



■f 

J — ( 



dtTr 



5L_d_ /SL 
SA It \Ja 



SA 



(10a) 



Assuming now that SA vanishes rapidly enough at t = ±oo so that all the Sqr vanish there, and that 
the generalized Euler-Lagrange equations are satisfied for all times, then both SS and the surface terms 
in Eq. (10a) vanish, and independence of SA{t) at different times implies 



= Tr 



S^^_d_ /SL 
SA dt \SA 



SA 



(10b) 



Suppose now that L is left invariant under the variations Sqr parameterized by SA. Then SL/SA = 0, 
and Eq. (10b) simplifies to 

, (10c) 



TV SA 



dt 
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with Qa the "charge" defined by 



Qa = 5 . (lOd) 
5A 



If L is invariant for arbitrary time-independent anti-self-adjoint operators 6A, then Eq. (10c) implies 
the operator statement 

n /in \ 

— = , (lOe) 

with Qa anti-self-adjoint. (Operator gauge transformations obey this condition trivially, with Qa = 0, 
since L is invariant for arbitrary time-dependent anti-self-adjoint 5A.) On the other hand, in the case 
of Poincare transformations, 6 A is a c-number describing an infinitesimal translation or proper Lorentz 
transformation of the coordinates, and Eq. (10c) then only implies the total trace relation 

^ TV Qa = , (lOf) 

of which Eq. (6d) is a particular example. 

An important special case is that in which 6A{t) parameterizes a linear transformation of the 

coordinates of the form 

Sqr = 6A{t) ^ Grsqs , (Ha) 

s 

with Grs independent of time and of the q's and g's. Then we have 

Sqr = 6A{t) Grsqs + 5A{t) ^ GrsQs , (Hb) 

S S 

which together with Eq. (3d) imphes that 

^G = 7T = ±E GrsQs^, (11c) 
6A 6qr 

rs 

again with the +(— ) sign chosen according as whether is of bosonic (fermionic) type. Substituting 
Eq. (5a) into Eq. (11c) and taking the trace, we thus get 

TrQG = Tr(^ PrGrsqs^ • (lid) 
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Let now Qh be the charge associated with a second hnear transformation in which Grs is replaced by 
Hrs', then for the generalized Poisson bracket of Tr Qg with Tr Qh, we find 



Hence if a set of matrices G,H, . . . used to generate Hnear transformations of the (/'s obeys a Lie algebra, 
then the corresponding functional Tr Qc TtQh, ■ ■ ■ obey the same Lie algebra under the generalized 
Poisson bracket. This shows, for example, that in a Poincare invariant theory defined by a total trace 
action, the total trace functional defining the Poincare generators will obey the Poincare algebra under 
the bracket operation of Eq. (6a). 
(4) When an operator gauge invariance is present, the problem of identifying physical observables becomes 
more subtle. Since observables should correspond to invariant geometric features of the quantum dy- 
namics, only operator gauge invariant quantities can be observables. Thus, the expectation values of 
Eq. (7b), which are not operator gauge invariant for general states |5), |/), are in general not observables. 
One way to form observable quantities is to construct total trace functional, similar to L but involving 
higher degree polynomials in the operators Qr, which are operator gauge invariant. Clearly, an infinite 
number of such observables can be constructed. A second way to form observable quantities is to focus 
on particular operators Ors transforming under operator gauge transformations as 



and on the co-transforming bases of states \br ), \bs ), \fr ), \ fs ),n = 1,2,..., which transform as 



{TrQG,TrQH} = Tr ^ {±){prGrs{±)Hstqt - PrHrs{±)Gstqt) 




(lie) 



(12a) 



(12b) 



Then the special class of matrix elements 



{bi^^\Ors\b<r^) , 



(/^"^laai/i")) 



(12c) 
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are operator gauge invariant, and hence observables. 
(5) As we have seen, in theories constructed from a total trace Lagrangian, the Hamiltonian dynamics 
is governed by the operator variational equations of Eq. (5f), which are generated by the total trace 
Hamiltonian H. There appears to be no special reason for this dynamics to be unitary, that is, in 
general there is no reason to expect that there should be a unitary time evolution operator U{t, 0), such 
that for all qr,Pr and all times t, the dynamics of Eq. (5f) is equivalent to 

q^{t) = U\t, 0)qr{0)U{t, 0) , Pr{t) = U\t, Q)pr{Q)U{t, 0) . (13a) 

An equivalent statement in infinitesimal form is that in general there is no reason to expect that there 
should be an anti-self-adjoint operator Hamiltonian H{t), such that for all qr,Pr and all times t, the 
dynamics of Eq. (5f) is equivalent to 

qr = [H{t),qr\ , Pr = [H{t),Pr] . (13b) 

It follows from these statements that total trace Lagrangian dynamics is, potentially, an even more 
general form of quantum mechanics than standard, operator Hamiltonian-based, quantum mechanics. 

Suppose, however, that the dynamics is such that a U{t, 0) and an H{t) obeying Eqs. (13a,b) do 
exist, for a theory with operator gauge invariance. Then we can make an operator gauge transformation 

qr{t) ^ U{t, Q)qr{t)U\t, 0) = qrs{t) , Pr{t) ^ U{t, 0)pr{t)U'l{t, 0) = Prs{t) , (13c) 

with a corresponding transformation for co-transforming states |6), |/), 

\b)^U{t,0)\b)^\bs{t)) , \f)^Uit,0)\f)^\fs{t)) , (13d) 

thus defining "Schrodingcr picture" operators qrs{t),Prs{t) and states \bs{t)) and \fs{t))- By construc- 
tion, the operators qrs{t) and Prs{t) are time independent, 

qrs{t) = qrs{0) , Prs{t) = Prs{0) , (13e) 

while the states \bs{t)) and \fs{t)) obey the Schrodinger time-development equation 

I \bs{t)) = -m\bsit)) , I \fs{t)) = -mmt)) . m 

14 



From this perspective, Schrodinger picture quantum mechanics appears as a rather special case of the 
more general quantum dynamics described by operator gauge invariant total trace Lagrangians. 

It clearly is of great importance to determine under what circumstances the operator time devel- 
opment equations of Eq. (5f) are equivalent to a unitary evolution as in Eqs. (13a, b). It seems likely 
that this equivalence always holds in complex quantum mechanics, since there the standard canonical 
quantization rules give a constructive procedure for going from the Lagrangian L, interpreted now as a 
classical Lagrangian, to an operator Hamiltonian obeying Eq. (13b). In the case of quaternionic quan- 
tum mechanics the situation is far from clear, and we leave the problem of determining whether and 
when Eq. (5f) is equivalent to Eqs. (13a,b) as an important open question. 
(6) As formulated up to this point, total trace quantum dynamics applies for arbitrary operator properties 
of the dynamical variables {5r}>{Pr} at some initial time (say t = 0), from which the Hamiltonian 
equations of motion can be integrated forwards to t > (and backwards to t < as well). Suppose 
now that for some subset of the variables {qr}, the Lagrangian L is independent of the time derivatives 
{cir}. Then the corresponding canonical momenta {pr} vanish identically. 



We are now dealing with a constrained Hamiltonian system, for which an operator generahzation of 
the standard Dirac treatment of constrained systems will be needed. This may involve some subtleties, 
since the general Dirac prescription employs second variational derivatives of the Lagrangian, which, as 
discussed above, are not well-defined in the operator case. The situation in which some of the canonical 
momenta vanish identically is of course not the most general form of a constrained system, but it is 
precisely what occurs when a gauge invariance is present. We assume, in analogy with the standard 
Yang-Mills case, that the correct procedure for operator gauge invariant systems will be to adjoin to the 
operator constraints of Eqs. (14a,b) an equal number of operator gauge-fixing conditions, which break 




(14a) 



and the Euler-Lagrange equations for these variables degenerate to the constraints 



5L 

Sqr 



(14b) 
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the operator gauge invariance. We conjecture that the constraints of Eqs. (14a, b), together with the 
operator gauge-fixing conditions, provide the minimum specification of operator properties of the {qr} 
and {pr} which are needed for a consistent theory. We will see, however, that there are examples in 
which it is possible to add further constraints beyond this minimum and still preserve consistency with 
the operator equations of motion. Finding a suitable operator generalization of the Dirac procedure, at 
least in the case of operator gauge invariant systems, and determining the precise conditions needed for 
operator specification, are again important open problems. 

(7) Finally, some historical remarks. The concepts of operator-valued gauge transformations and a total 
trace action were introduced, without the (—1)^ factor and real part of Eq. (lb), by Adler [3,4,5] in 
the context of a theory termed "algebraic chromo dynamics." The (—1)^ factor was used by Witten 
[2] to define a topological index, which in our notation is simply Trl. The necessity for including 
the real part in the definition of a trace for quaternionic Hilbert space was pointed out by Finkelstein, 
Jauch, and Speiser [6] . A general suggestion of operator-valued gauge transformations was also made by 
Mackey [7], with an implementation [8] which formulates the equivalence class, under unitary operator 
transformations, of complex Galilean invariant Hamiltonians. Since what we have done above is to set 
up a dynamics on a manifold with non-commuting coordinates {qr}, the discussion of this and the 
subsequent two sections appears to be related to the non-commutative geometry program of A. Connes 
[9,10]. There are close analogies between the identification of observables in operator gauge invariant 
theories, and the identification of observables in general relativity and in conventional Yang-Mills gauge 
theories. Operators transforming as in Eq. (12a) are analogs of bitensor quantities in classical general 
relativity and of path ordered integrals in gauge theories. The suggestion that the generalization from 
special to general relativity should have an operator analog in the generalization from complex to 
quaternionic quantum mechanics, was first made by Finkelstein, Jauch, Schiminovich, and Speiser [11]. 
Because the analysis of this section has dealt with the general case, it has of necessity been rather 

abstract. Concrete illustrations of operator gauge invariant systems, in complex and quaternionic quantum 

mechanics and quantum field theory, are given in the next two sections. 
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3. Operator gauge invariant total trace Lagrangian formulation 
of complex quantum mechanics 

We proceed in this section to illustrate the general formalism which we have just set up, in the familiar 
context of complex quantum mechanics and quantum field theory. In complex quantum mechanics, the 
left-acting operator / defined by 

I = Y,\n)i{n\ (15) 

n 

commutes with all operators. (This permits the identification of / with il, with 1 the unit operator and i 
the right-acting imaginary unit, a notation which is standard in the complex quantum mechanics literature 
but which we shall not follow here. A discussion of right-acting versus left-acting algebras is given in 
Sec. 4 below.) Hence the complex specialization of operator gauge invariance consists in assuming that the 
operators Ur, U'^ of Eqs. (8a) and (9a,b) commute with 7, which permits the inclusion of explicit factors of 
/ in the construction of operator gauge invariant total trace Lagrangians. 

As our first example, we consider a single self-adjoint bosonic coordinate q{t) obeying Galilean invariant 
dynamics. The conventional Lagrangian for this model (with the mass taken as unity for convenience) is 

L'',^\e + \{q,A{q)}+-V{q) , (16a) 

with A{q) and V{q) self-adjoint functions of q and { , }+ the conventional anti-commutator; this Lagrangian 
clearly does not have any simple transformation properties under the operator transformation 

q-^UqU\ UW^U^U=1. (16b) 

To achieve covariance, we follow the standard procedure of replacing the ordinary time derivative do = d/dt 
by a covariant derivative Dq, defined as 

Doq=^ + [Bo,q\, (16c) 

with Bq an anti-self-adjoint operator gauge potential. Under the transformation of Eq. (16b), Bq is taken 
to transform as 

Bo ^ UBoU'^ _ [/t = c/BoC/t + [/ [/t (i6d) 
at at 
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as a consequence of which Doq transforms as 



= u(^^ + [Bo, q]^ = U{boq)U^ . (16e) 



Hence if we redefine the Lagrangian of Eq. (16a) as 



1 .2 1 



L, = - {Doqy + - {Doq, A{q)}+ - V{q) , (17a) 
then under operator gauge transformations Lg transforms covariantly, 

Lg ^ [/Lgf/t , (17b) 
and the corresponding total trace Lagrangian and action defined by 

/oo 
dtLg , (17c) 
-oo 



Lg Tr Lg , Sg 



are invariant. 

We have now achieved operator gauge invariance at the price of introducing an extra dynamical variable 
-Bo- We must next investigate the structure of possible Lagrangians to govern the dynamics of Bq- The 
usual Yang-Mills type field strength Fqo vanishes identically, 

Foo = doBo - doBo + [Bq, Bo]=0, (18a) 

and so cannot be used to construct a Lagrangian for Bq. However, there is one additional Lagrangian which 
can be formed from Bo, 

LBo=Tr{IBo), (18b) 

with the inclusion of a factor of / necessitated by the fact that Bq is anti-self-adjoint. To see that Eq. (18b) 
defines a satisfactory Lagrangian, we note that it suffices to check its behavior under infinitesimal operator 
gauge transformations of the form 

U = 1 + 5A, 6A = -Mt , (18c) 
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under which the first order variation of Bq is 



5Bo = [5A, Bo]-^ = -DoSA . (18d) 



From Eqs. (18b,d) we find 



5Lb, = Tr {I6B0) = TV ([Bq, I]5K) - ^ ^ (^'^A) = Tr (7M) , (18e) 
which is a time derivative, and as a consequence the total trace action S_Bo defined by 

/oo 
rfiLso (19a) 
-00 

is invariant under Eq. (18c) when 5A vanishes &i t = ±00 [or more generally, when (5A(oo) = 5A(— 00)], 

5Sb, = . (19b) 

This argument of course does not imply that Sb^ is invariant under the global transformation of Eq. (16d), 
which can produce changes in the topological sector; rather, what we have shown is that S Bo defines a form 
of topological action, which is constant within each distinct topological sector defined under operator valued 
gauge transformation. 

An alternative argument for the invariance of Sbo proceeds directly from Eq. (16d), which implies that 

Tr (/Bo) ^ TV 7 (UB<,U^ ~ ^ = ^ (^^0) ~ ^ ^0 ' ^^^""^ 

Writing U = e^^^\ with A = — A^, and using the operator identity (which can be verified by power series 
expansion and term-by-term integration) 

^ e^- dse«^^e(i-«)^, (19d) 



dt 7n dt 

we have 



Hence Sbq is invariant whenever Tr [/A(oo)] = Tr[/A(— 00)]. This argument, as well as the infinitesimal 
one given above, extends immediately to the case in which Bq is replaced by a space-time component of a 
four vector gauge potential and d/dt is replaced by a space-time derivative. 
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(20b) 



Let us now examine the dynamics following from a general linear combination of the actions and 
Sbq- Forming the total Lagrangian L and action S, 

L = Lg — AqLbo , S = Sg — AoSsg , (20a) 

with Ao a constant, and taking general operator variations, we get (with repeated use of the cyclic property 
of the trace) 

6L = Tr {[^og + A{q)]S{Doq) + iDoq)dA{q) - SV{q) - XoISBq} 

= Tr {[Doq + A{q)]iSq + [SBo, q] + [Bo, Sq]) + F{q, Doq)Sq - XqISBo] 
= Tr { [^05 + A{q)]dq + {[Doq + A{q), Bo] + F{q, Doq)) 6q 
+ ([q,Doq + Aiq)]-Xol)6Bo} , 
where we have defined a generalized force term F{q, Doq) by 

Tr [{boq)5A{q) - 5V{q)] = Tr [F{q, Doq)Sq] . (20c) 



Hence 

— =Doq + A{q), — = [Doq + A{q),Bo]+F{q,boq), (20d) 

^=0, ^ = [q,boq + Aiq)]-XoI , (20e) 
5Bo oBo 

and so the Euler-Lagrange equations following from (5S = consist of a dynamical equation for q{t), 

bo[boq + A{q)] = ^ [Doq + A{q)] + [Bo, Doq + A{q)] = F{q, Doq) , (20f) 
together with a constraint 

[q,boq + A{q)]=XoI , (20g) 

both of which are covariant under the operator gauge transformations of Eqs. (16b-e). When we rewrite the 
dynamics in total trace Hamiltonian form, we identify the canonical momentum conjugate to q as 

P=Y^ = boq + A{q) , (21a) 
dq 
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and so the constraint of Eq. (20g) reads 

[q,p]=XoI, (21b) 

and is just the standard canonical commutator with the identification Aq = h{= 1 in microscopic units). 
Note that because of the real part in the definition of Tr, we have Tr/ = 0, and so taking Tr of the left- 
and right-hand sides of Eq. (21b) gives a consistent equation for non-zero Aq, 

= Tr [g, p] = AoTr 7 = AoO . (21c) 

Carrying out the Legendre transform of Eq. (5b) we get 

= Tr Uooq + A{q)]{Doq - [Bq, q]) - \ {D^qf - {Doq)A{q) + V{q) + Aq/Bo 



= 1^ |i {Doqf + V{q) + {XoI-[q,Doq + A{q)])Bo^ 



= TV |i [p - A{q)]^ + V{q) + {Xol - [q,p])Bo^ , (21d) 
which in the "Hamiltonian gauge" 

Bo = (21e) 



simplifies to 



H = Tr{i[p-^(9)]^ + y(g)}. (22a) 



Taking the operator variation of Eq. (22a), and recalling the definition of Eq. (20c), we get 

<5H = Tr {[p - A{q)][Sp - SA{q)] + SV{q)} = Tr {[p - A{q)]Sp - F{q,p- A{q))Sq} . (22b) 

Hence the total trace Hamiltonian equations of motion are 

q=-^=p-A{q), p=——=F{q,p-A{q)), (22c) 

which agree with Eqs. (16c), (20f) and (21a) when these are specialized to the gauge Bq = 0. 

We see, then, that the Heisenberg picture equations of motion and the canonical commutation relation 
for a Galilean particle both emerge from the operator gauge invariant, total trace Lagrangian formalism. The 
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derivations just given are expressed in operator terms throughout; at no point did we introduce a classical 
Lagrangian and its "quantization." In complex quantum mechanics, the conventional canonical quantization 
route is of course still valid, and implies that there is an operator Hamiltonian H given by the Weyl ordering 
of the corresponding classical Hamiltonian. For the model under study, we have 

H = {\\p-A{q)f + V{q)}w, (23a) 

where the Weyl ordering subscript W implies symmetrization of p with respect to the factors of q in each 
term of A{q), e.g., 

= {(tp}w = ^-r (O + + 9"- W + • • • + gV""' + w""' + w") ■ (23b) 

n + 1 

With {q,p] = XqI, we then find that 

±-[q,H]=p-A{q), -L. h] = -F{q,p - A{q)) , (23c) 

since, as may be verified by some algebra, the Weyl ordering of H leads to the same factor ordering in the 
force term F as is obtained, via the cyclic property of the trace, from the operator variational definition of 
F in Eq. (20c). Hence the equations of motion of Eq. (22c) are equivalent to 

q = Xo'[IH, q], p = X^'ilH, p] , (23d) 

and so in the terminology of the preceding section, the dynamics is unitary. This permits us to transform 
from the Heisenbcrg picture to the Schrodingcr picture, in which the operators are time independent and 
the co-transforming states carry the quantum dynamics. 

In the example just given, all operators are bosonic, and so the (—1)^ factor in the definition of Tr does 
not come into play. As our second example, we consider a single non-interacting fermion degree of freedom 
with mass m, described by the conventional Lagrangian 

= ^ (^t^ _ ^t^) _ TO^t^ . (24a) 
The Lagrangian of Eq. (24a) again does not transform simply under the operator transformation 

V' ^ UipU^ , UU^ = U^U = 1 , (24b) 
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but as in the bosonic example, we can achieve covariance by replacing the time derivative by the covariant 
derivative Dq, 

Do^=^ + [Bo, i>] , {Do^)'< = Do^'< = ^ + [Bo, V'^] • (24c) 
With Bo transforming as in Eq. (16d), Dotp transforms as 

Doi^ U{Doi})U^ , Doi^^ ^ U{Do^^)U^ , (24d) 

and so the redefined Lagrangian 

Lv, = ^ ii'^Doi) - (^oV'^)V'] - mV'V (24e) 

transforms covariantly, 

Lv- ^ UL^U^ . (24f) 
The corresponding total trace Lagrangian and action , 

/oo 
dt , (24g) 
-oo 



/"OO 

-oo 



are then invariant. 

Proceeding as in our first example, we examine the dynamics following from a general linear combination 
of the actions and Sbq- Writing 

L = — AqLbo , S = — AqSbq , (25a) 

and taking general operator variations of ijj and Bq, with Stp^ = (5^))^ we get 

(5L = Tr 1^ 5jI>''DojI> + V''^(5(£>oV') - ^(£'o^^)V' - {Doi'^)Sil> 
-to((5^/'V + V-t^V) - Ao/^So} 
= Tr (2{Si;)^{IDoiJ - mV^) - /[{^, + Ao](5Bo) + ^ Tr {1^;^^;) . (25b) 

We have here explicitly used properties (ii) and (iii) of Tr discussed in the preceding section, with the 
appearance of the anti-commutator {^, •il>^}+ a direct result of the effect of the (—1)^ factor on the reordering 
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of fermion factors inside Tr. The time derivative term in Eq. (25b) makes no contribution to the action S, 
and so the Euler-Lagrange equations following from (5S = can be read off directly from Eq. (25b), giving 
the dynamical equation for ip 

IDoip = mip , (25c) 



together with the constraint 



{V',V^}+ = -Ao 



(25d) 



Equation (25d) is just the standard canonical anti-commutator for a fermion degree of freedom when we 
identify — Aq = 7i (= 1 in microscopic units). We note that this determination of Aq has the opposite sign 
from that found above in the bosonic example; we shall say more about this shortly. Finally, we remark that 
because the Hilbert space for a single fermion degree of freedom consists of a zero fermion state |0) G 
and a one fermion state |1) G V^, the effect of the (—1)^ in the definition of Tr is to give 



Trl = l-1 = , 



(25e) 



and so taking Tr of the left- and right-hand sides of Eq. (25d) gives a consistent equation for non-zero Aq, 



= Tr (^V^ - ^V'^) = Tr {V', = -AoTr 1 = -AqO . 



(25f) 



From the time derivative term in Eq. (25b), we identify the momentum canonical to ip as 



dtp 



Since by property (ii) of Tr we have 



Tr(/V'^V') = Tr 



the total trace Hamiltonian becomes 



- tp^ 



Tr 



- {11)^ - V^V) 



(26a) 



(26b) 



H = Tr (/V'^V') - L = Tr - {tpH - tpH) 

= Tr (/[{^, + Ao] Bo + mVV) , 
24 



-L 



(26c) 



which simpUfies in the Hamiltonian gauge Bq = to 

H = Tr (m^/^V) = Tr {-Imp^ip) . (26d) 

Hence taking operator variations, we get 

(5H = Tr {—Im dp^^ — Imp^S^) = Tr [Im {ip6p^ — p^5tp)] , (26e) 

and so Eq. (5f) becomes 

lip'' = pw, = — — = Imp^ = -rml)^ , V' = — = -Imip , (26f) 
Sip ^ Sp^ 

in agreement with Eq. (25c) and its adjoint when Bq = 0. As in our bosonic example, in the fermionic case 
there is also an operator Hamiltonian 

H = m^j'^ip (26g) 

which, with — — Aq, obeys 

^ [i;,H] = -Imtp = ip , [ip\ H\ ^ Imtp^ ^ . (26h) 



-/Ao ^ " " ' -/Ao 

Hence the dynamics of Eq. (26f) is unitary, again permitting an operator gauge transformation from the 
Heisenberg picture to the Schrodinger picture. 

We see, then, that the standard quantum mechanics of a single bosonic and a single fermionic degree 
of freedom, including the canonical commutator and anti-commutator, follows from the operator gauge 
invariant total trace Lagrangian formalism. Let us now examine what happens when more than one degree 
of freedom is present. Focusing only on the structure of the kinetic terms, for a set {qr}, {V's} of bosonic 
and fermionic variables, we have the conventional Lagrangian 

R S 
r— 1 s—1 

There is now more than one way to extend Eq. (27a) into an operator gauge invariant Lagrangian, depending 
on whether we require invariance when all variables are subject to the same operator unitary transformation, 
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or are subject to independent operator unitary transformations. When all variables are subject to the same 
transformation, 

qr UqrU^ , ^ Uij.U'' , UU^ = U'^U = 1 , (27b) 

we achieve covariance of the Lagrangian by replacing L'^^^y ^^^^ by i-^g],} {,/;,}' with 

^{luM = 11 I (^oqrY + E ^ [V'l^oV'. - {Doi'D^s] , (27c) 

r=l 5=1 

where Dq acts on the bosonic and fermionic degrees of freedom as in Eqs. (16c) and (24c). The total trace 
Lagrangian and action 

/oo 
-OO 

are then invariant under Eq. (27b), and the constraint equation arising from the variation with respect to 
Bo of 

^ = ^MA^s] ~ ^o^^o > (27e) 
with Sbo defined as in Eqs. (18b) and (19a), is [cf Eqs. (20b) and (25b)] 

R S 

[qr, Doqr] -lY, {i^s,4}+ - Ao/ = . (27f) 

We see that the individual canonical commutators and anti-commutators are not determined, but only the 
linear combination given by the sum of bosonic commutators minus the sum of fermionic anti-commutators. 
The constraint of Eq. (27f), as well as the dynamical equations of motion for {qr} and {tps}, are consistent 
with the imposition of the canonical relations 

[qr, Doqr] = I , {i's, V'|}+ = 1 , (27g) 

for each r and s, but do not require these. In the language of the conventional theory of constrained systems, 
Eqs. (27g) are invariant relations which are compatible with the Lagrangian constraints and equations of 
motion. We note the interesting fact that when the numbers of bosonic and fermionic degrees of freedom 
are equal, as is the case for a supersymmetric theory, and the canonical relations of Eq. (27g) are imposed, 
then the constraint of Eq. (27f) is satisfied with Aq = 0. 
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An alternative possibility is to require invariance of the total trace Lagrangian under independent 
operator transformations of the canonical variables, 



UrUl = UlUr = UsUl = UlUs = 1 . 



(28a) 



To achieve this, we introduce an independent covariant derivative for each canonical variable. 



(28b) 



with Bqj., r = 1,. . . ,R and Bqs, s = 1,. . . ,S independent anti-self-adjoint operator gauge potentials. We 
now replace the Lagrangian L'^g^j ^ by 



R S 



(28c) 



the individual terms of which transform covariantly (but now independently) under Eq. (28a), 



r=l 



ul. 



(28d) 



Equation (28d) has the general form of Eq. (8d), and correspondingly, the total trace Lagrangian and action 
defined by 

/oo 
"^^^llu^s} (28e) 
-oo 

are invariant under Eq. (28a). For each gauge potential Bor,s, we can add an action term S^or.s formed as 
in Eqs. (18b) and (19a), giving for the overall total trace action 



(28f) 



r=l 



Varying with respect to each -Bor,s, we get the independent constraints 



[qr, Dorqr] - Aor-f = , r = 1, . . . , R 

{i;s,i>l} + Xos = , s = l,...,S., 
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(28g) 



which are closer in structure to the usual canonical commutators than the single constraint of Eq. (27f). 
Of course, when we impose the requirement of invariance under independent operator transformations as in 
Eq. (28a), the interaction terms in the total trace Lagrangian are much more tightly restricted in form than 
when we impose only the global invariance of Eq. (27b) . In field theory applications, the indices r, s are 
typically composite indices, indicating both the spatial coordinate value x and the particular field component 
at X. In this case, the constraints associated with Bq have a structure intermediate in form between those of 
Eqs. (27f) and (28g): the constraints at different values of x are independent, but at each x consist of a sum 
of contributions from the various bosonic and fermionic field components present in the theory, evaluated at 
that value of x. 

Various further examples of operator gauge invariant quantum mechanical and quantum field systems 
in complex quantum mechanics, including models with independent left and right gaugings as in Eq. (9b), 
and operator gauge invariant extensions of Yang-Mills theory, are given in [1]. In fact, we show there that all 
of the basic field theory building blocks of the standard model can be imbedded in operator gauge invariant 
theories; this raises the question of studying operator gauge invariant extensions of the full standard model, to 
see if useful insights (such as restrictions on the parameters, or new calculational methods) can be obtained. 

4. Operator gauge invariant quaternionic field theories 

We turn now to the construction of operator gauge invariant theories in quaternionic quantum mechanics, 
which is distinguished from the complex case discussed in the preceding section by the non-existence of a 
left-acting I which commutes with all operators. For what follows only a few facts about quaternionic 
Hilbert space are needed, which we now proceed to state. Quaternionic Hilbert space Vh is closed under 
taking linear combinations of ket vectors when right multiplied by quaternion scalars. A quaternion scalar 
is a number q of the form 

q = qo+iqi+ jq2 + kqs , go,i,2,3 real , 

= = = —1 , ij = —ji = k , jk = —kj = i , ki = —ik = j , (29a) 

from which one finds that quaternion multiplication, while non-commutative [i.e., qiq2 7^ q2qi in general], is 
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associative [i.e., 91(9293) = (9192)93]- The conjugate 9 and modulus \q\ of a quaternion are defined by 

q = qo-iqi- jq2 - kq^ , 

|9| = (9-9)'/' = (99-)'/' = {qI + qI + qI + , (29b) 

and by Eq. (29a), the conjugate of a product is the product of the conjugates in reverse order, 

9i92 = 9291 • (29c) 

Defining the adjoint in Vh as the quaternion conjugate of the matrix (or operator) transpose, the adjoint 
of a ket vector is a bra vector, the adjoint of a scalar is its conjugate, and from Eq. (29c) one finds the 
usual rule that the adjoint of any product of multiple factors is the product of the adjoints of the factors in 
reverse order. The inner product {f\g) is quaternion-valued, and so any two inner products {f\g) and {f'\g') 
do not in general commute with one another in quaternionic Hilbert space (whereas in complex quantum 
mechanics they do). In other words, in quaternionic Hilbert space not only do operators not commute with 
one another, but their matrix elements as well are non-commutative. 

For any complete orthonormal set of states {|n)}, one defines [12,13] the anti-self-adjoint left-acting 
algebra I,J,K by 

{I,J,K) = J2 |n)(i,i,fc)(n| = -(/t,jt,i^t) , (29d) 

n 

which obeys an isomorphic image of the algebra of right-acting quaternion scalars, 

J2^A'2 = -1 = -^ |n)(n| , 

n 

IJ = -JI = K , JK = -KJ = 1, KI = -IK = J . (29e) 

For an arbitrary quaternion operator O which acts on kets from the left (all quaternion quantum fields 
discussed below fit this description), we can introduce [12,13] a set of "formally real" components Oa, A = 
0,1, 2, 3 defined by 

00 = 1 (O - lOI - JO J - KOK) , 

01 = - J {lO + 01- JOK + KOJ) , 
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02 = --{J0 + 0J- KOI + lOK) , 

Og = - 1 (KO + OK- lOJ + JOI) , (29f ) 

so-called because they obey 

O = 0o + I0i + J02 + K0z, [Oa, {I,J,K)]=0 , A = 0,l,2,3, (29g) 

but they do not in general commute with one another. Finally, we note that — ^ "^'"^ ^ ^'^^ ^ 

anti-self-adjoint generators of a one-dimensional quaternionic irreducible representation [6,14] of the group 
SU{2); since the smallest non-trivial complex irreducible representation of SU{2) is the two-dimensional 
spinor representation, this means that any quaternionic quantum field theory based on gauging the one- 
dimensional quaternionic representation of SU{2) is irreducibly quaternionic, and cannot be reduced to a 
complex quantum field theory. All of the field equations discussed in this section have this character. 

With these facts in mind, we are ready to begin the construction of operator gauge invariant quaternionic 
field models. We begin with the field theory of a quaternionic scalar field (p, which is not restricted to be 
self-adjoint (or anti-self-adjoint), and which is subjected to independent left and right local gaugings,* 

^ C/</>[/'1' , UU^ = U^U = U'U'^ = U"'U' = 1 . (30a) 
Introducing anti-self-adjoint gauge potentials Bfj_,B'^ which transform as 

^ UB^^U^ - {d^U)U^ , ^ U'B'^U'^ - {d^U')U'^ , (30b) 

and the covariant derivative and field strengths 

D^<i) = d^(i> + B^<i>-(j>B'^ , 
F^v = d^Bi, — d^Bn + [B^, Bi] , 

Fl^ = d^Bl-d.B'^ + [B'^,Bl], (30c) 
* An historical discussion of various gaugings which have been proposed for quaternionic fields, and 
references, is given in Sec. 11.2 of [1]. 

Our metric convention is goo = — 1, 9ii= 922 = 933 = 1- 
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which transform as 

D^<f>^U{D^<f>)U'U F^^^UF^^U^, Fl^^U'Fl^U'^ , (30d) 
the operator gauge invariant total trace Lagrangian density for the scalar field is 

C = jlLif, + Cb + C,B' , 



The total trace Lagrangian L and action S are formed from C by the usual recipe 

1, = j d^xC , S = j dtL . 



(30e) 



(30f) 



Note that an action term analogous to Sbq of Eq. (19a) is now not admissible, because the left-acting I 
needed to construct this term breaks the operator gauge invariance. When we vary S, through SF/^^, and 
SF'^^, we encounter the covariant derivatives and D'^, 

SF^, = D^dB, - b,5B^ , 5Fl, = %SBl - D'JB'^ , (31a) 
and in integrating by parts we use the intertwining identities [1] 

d^Tr (prj^) = Tr [{D^pW + - 

d^Tr ip^rj) = Tr [iD^p)^ + p^D^rj] , (31b) 

with p and r] either both bosonic or both fermionic in type. Omitting further computational details, we get 
the operator equations of motion 

D^Di'cf) - (m^ + X(j)(f)^)(f) = , 

Di'F,^ = G^X , J. = \ [0(1^.0)^ - (i?.0)</>1 , 

D'^Fl^ = G'^X , X = l [<I>'^D.<I> - (Z).</.)t0] , (32) 
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in which the u = components of the gauge field equations are constraints rather than dynamical equations. 
Equations (30e) and (32) can be specialized to less general gaugings of 0. For example, if we take </> to be 
gauged under 

(j) U(j)U'i , (33a) 
which is the most general allowed gauging when (f) is self-adjoint, the appropriate Lagrangian density is 

C = C^ + Cb, = IV |i [-(^^</.)tD''</> - m20t<^] _ ^ (<^t<^)2| ^ (33^^ 

and the corresponding equations of motion are 

r>^Z)^0 - (m^ + X(l)(j)'')<f) = , 
D'"F^^^^ = \g'' [ci){D.(l))^ - {Du<f>)(p'' + 4>'^b.(j> - (^.<^) V] • (33c) 
We turn next to the case of quaternionic fermion fields, working always in the Majorana representation 
for the Dirac matrices, starting again with the most general gauging, in which there are two fermions V'(i) , V'(2) 
transforming under independent left and right local gauge transformations as 

V(i) ^ f/V(i)«7'^ , ^(2) ^ t^^(2)f/'^ . (34a) 
The operator gauge invariant total trace Lagrangian density is 

C = „ +Cb + Cb' , (34b) 
with Cb and Cb' as in Eq. (30e), and with C^^^ given by 



^ V';2)7°7^^,.V'(i) + (i5;.V(i))^7°7^V'(2) - V;i)7°7'^^mV'(2) - (l?^V(2))W°7^V'(i) 
+m 



(v.;2)Z7V(i)-V'fi)nV(2))} • (34c) 
Varying the action S [still related to C by Eq. (30f)] and recalling that in the Majorana representation 7^7'* 
and i'y^ are, respectively, real symmetric and real skew-symmetric matrices, we get the operator equations 
of motion 

(707''£>^ + mi7°)^(i) = , {l^'l'^D^ + mi7°)V(2) = , 
b^'F,^ = G^J, , X = V'fi)7j7°^V'f2'i - i^i2)lll°^4i) > 

^'"F^^ = -{G'fX , X = ^^;i)7°7.V'(2) - V'f2)7°7.^^(i) , (35) 
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with T indicating Dirac index (but not operator) transposition. Again, the y = Q components of the gauge 
field equations are constraints. As in the boson case, we can readily specialize the two fermion model to the 
less general gauging 

^(i) ^ C/^(i)C/^ , ^(2) ^ U^(2)U^ , (36a) 
for which the appropriate Lagrangian density is 



[I [V'f2)7°7''£>MV'(i) + (I);.V(i))^7VV'(2) - Vfi)7V^^V'(2) - (^^V(2))^7°7''V'(i) 
+m (^/;J2) i7°V'(i) - V'(i) «7°V'(2)) [ , (36b) 



and the corresponding equations of motion are 

(7°7''D^ + mn°)i/'(i) = , (7°7''-D;. + mi7°)V(2) = , 



D^^F,^ = [Vfi)7j7°^^S - ^(2)lul°^^li) - V';i)7°7.V'(2) + ^;2)7°7.^(i) • 
We can also form fermionic models with a single fermion field, in which either the left or the right gauge 
invariance is restricted to be a complex gauge invariance. With the right gauge invariance restricted, we get 

= Tr 1^ [V^ 7°7''£»^^ - (£>^V')H°7''V'] + /'m^ti^o^l ^ (37^) 

with I' a space-time independent left algebra operator, and with U' and B'^ restricted to be C(l,/'). 
Similarly, with the left gauge invariance restricted, we have 

£; = TV |i [iP^Ij^jf^D^^ - (i?^V)W°7''/V'] + mV+/n>| , (37b) 

with / a space-time independent left algebra operator, and with U and i?^ restricted to be C(l,/). When 
is so restricted, an action term Ss^ analogous to Eq. (19a) can be included in the total action, and 

similarly for B'^ in the case of Eq. (37a). 

Before proceeding to the total trace Hamiltonian form of the dynamics, we discuss a number of issues 

which can be addressed directly from the total trace Lagrangian and the equations of motion. 
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(36c) 



(1) We begin by contrasting the quaternionic gauge field structure with that of a conventional Yang-Mills 
gauge field. Let 1,Ea, A = 1,2, 3 be a space-time-independent left algebra basis, and let us use 
Eqs. (29f,g) to expand the gauge potential and the corresponding field-strength F^f^y over this basis, 

3 3 

= + ^ Ba^Ea , F^i, = i^o/ii/ + ^ Fa^vEa , (38a) 

with the expansion coeflacients Bo/i, BAf^, Fofj^^,, FAfiu formally real, 

[Bo^, Ec] = [BA^,, Ec] = [Fo^., Ec] = [Fa^., Ec]=0. (38b) 

We recall, however, from the discussion beginning the section, that in general the expansion coefficients 
Bo(n, Bai^, ■ ■ ■ are still operators which do not commute with one another. Substituting Eqs. (38a) into 
the formula of Eq. (30c) which relates the gauge field strength to the gauge potential -B^, we find 
that the expansion coefficients Fa^u, ^ = 0, . . . , 3 are related to the Ba^, ^ = 0, . . . , 3 by 

3 

Fq^v = d^BQi, — OvBq^ + [Bqh, Bq^] - ^ [Bah, Bav] , 

(38c) 

FAfiv = Q^Bav — dvBAfji + [Ba^, -Boy] - [Bav, Bq^] + ^ SABciBsij,, Bcv}+ ■ 

B,C=1 

If Bah, ^ = 0, . . . , 3 all commute with one another, Eq. (38c) would reduce to U (1) and SU{2) conven- 
tional gauge field structures, 

Fon = dfiBdu — dvBdfj, , 

a 

Fa/xv = O^Bav — dvBAn + 2 ^ £abcBbiiBcv ■ (38d) 

B,C=l 

But in the general case with non-commuting formally real components B^^, Eqs. (38c) are not equivalent 
to Eqs. (38d). Equations (38c) represent only part of the complete system of equations following from 
the total trace Lagrangians of Eqs. (30e), (34b, c), etc. It is straightforward to reexpress all of the 
remaining field equations in terms of formally real components with respect to the left-acting algebra 

1,Ea. 

Because the Ba^ are quaternionic operators, they can themselves be expanded over formally real 
components with respect to a second left-acting algebra which commutes with 1,Ea, 

3 

Baii = BqAh + ^ Bba^E^^^ , 

B=l 
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[Bbaij.,Ec] = [Bbah,E^^] = , 
[Eb,E^^^]=0, A,B,C = 0,1,2,3. 



(39a) 



This process can be continued to any order, giving in nth order a multi-quaternion expansion of the 



form (with Eq = Eq 



(1) 



E, 



(n) 


3 3 



1) 



^.= E E •••E B^_^,^,E^E^l..E^-^, 

A=0 Ai=0 A„=0 



(39b) 



with, for aU index values, vanishing commutators 



[Eb,eP] = [4^\4^^] = [Ba„.m^,Eb] = [Ba^...ma^.,E'^^>] = . 



Ml 



(39c) 



Note that the occurrence of left acting multi-quaternion algebras does not imply that probability am- 
plitudes belong to a non-division algebra. The algebra of right-acting scalars remains just 1, k, with 
the possibility of multiple commuting left-acting algebras arising because these are quaternionic Hilbert 
space operators, rather than scalars. For example, in as small as a four dimensional quaternionic Hilbert 



space there arc two miitiially commuting left algebra bases, 

,Ei=il, E2=jl, E3 = kl , 



E 



1 



(1) 



10 

10 

10 

1 
1 
0-10 
10 

-1 0, 



10 
1 
-10 
-10 



0-10 0' 
10 
1 
0-10. 



(39d) 



Properties of multi-quaternion algebras have been studied in a series of papers by Razon and Horwitz 
[15-17], and by various authors [18-21] who classified the U{2) case of the "color charge algebras" 
introduced by Adler [22] . (The color charge algebras correspond to multi-quaternion bases E^^^^ . . . E^^^ , 
in which the indices Ai,. . . ,An are contracted so as to leave a single free index A.) The latter calculations 
indicate that large Lie algebras are readily built up from multi-quaternion bases. 
(2) In the fermionic Lagrangians of Eqs. (34c) and (36b), the gauge bosons couple to vector currents 
only, and so there are no chiral anomalies and also no Witten [23] anomalies. Since the 75 matrix in 
Majorana representation is imaginary, attempting to split Eqs. (34c) and (36b) into chiral components 
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would break the operator gauge invariance. In other words, insisting on a maximal operator gauge 
invariance in quaternionic field theory excludes chiral fermions. 

At the same time, since 175 is real and anti-self-adjoint, when the mass m is zero the fermion total 
trace Lagrangians written above are invariant under 

V'(l,2) - , V;i,2) - ^l,2f~'''^ ' 

with f3 a real c-number which is independent of x. Thus, when m = 0, the fermion models all have a 
chiral symmetry. 

(3) Let us next address the issues of the spin-statistics connection and discrete symmetries in quaternionic 
field theory. On spin and statistics, we have little to say, beyond the fact that the conventional spin- 
statistics connection has been assumed in including the (—1)^ factor in the definition of Tr. We have 
not made a study of the possibility of an abnormal spin-statistics connection, but note that this would 
lead to the usual pathologies upon specializing a quaternionic field theory back to a complex one. 

To study discrete symmetries of the Lagrangians written down above, we use some standard Dirac 
matrix apparatus. We begin with parity (P), and readily find that under the substitutions 

Bo{x,x°) ^ Bo{-x,x°) , B'o{x,x°) ^ B'o{-x,x°) , 

Bi{x,x°) -Be{-x,x°) , B'i,{x,x°) -B[{-x,x^) , 

(j){x,x°) 'qp4>{-x,x°) , 

V'(i,2)(^,a;°) ^ Tl'pil^^(i,2){-x,x^) , ipix^x^) r]pij°'ii;i-x,x°) , (40a) 

with r]p,r]'p, and rjp arbitrary real c-number phases, all of the total trace Lagrangian densities of this 
section transform as 

C{x,x°) ^ C{-x,x°) , (40b) 

and the corresponding total trace actions are invariant. Turning next to time reversal (T), the relevant 
substitutions are now 

Bo{x, a;°) ^ -Bo{x, -a;°) , B'o{x, a;°) ^ -B'o{x, -x^) , 
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Be{x, x°) ^ Be{x, -x^), B[{x, :c°) ^ B'^{x, , 

V'(i)(x,a;°) ^ , 

ij)(2){x,x^) -77rAV'(2)(^, -a;°) , 

V'(x,a;'') ^ <^ /, ON r z'' > (41a) 

( ry^ JA^(x, for , ^ ^ 

with r]T,i]'T, Vt arbitrary real c-number phases, with A the Majorana representation time reversal ma- 
trix, which is real, and with J', J anti-commuting respectively with /', I. Under these substitutions all 
of the total trace Lagrangian densities of this section transform as 

C{x,x°)^C{x,-x°) , (41b) 

and the corresponding total trace actions are again invariant. 

Turning finally to charge conjugation (C), we consider the substitutions 

V'(i) v'ci^li) , V'(2) ^ -v'ci^l2) . 

V - vc^^^ , (42a) 

under which the covariant derivatives £'y:j0, D^'!/'(i,2)) and Di^ip transform as 

D^^ ^ VciD^cj,)^ , D;x^(i,2) ^ (+, -)r?^(D^V(i,2))^^ , D^^j ^ VciD^i^)^^ , (42b) 

with r]c,il'c^ ^^'^ v'c arbitrary real c-number phases and with T the Dirac index transpose. When we 
impose a condition of equality on the gauge field couplings, 

G = G' , (42c) 

all of the total trace Lagrangian densities of this section are invariant under the substitutions of Eq. (42a), 

C^C . (42d) 
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Since the gauge potentials and B'^ are interchanged by the substitutions of Eq. (42a), they are not 
C eigenstates. We see, then, that when the requirement of C invariance is imposed, the models with 
independent left and right gaugings are left with a single coupling constant G* 
(4) By using the fact that 175 is real and anti-self-adjoint, we can construct a model with self-adjoint 
Lagrangian using a single fermion field, without breaking the bi-unitary operator gauge invariance. 
The Lagrangian density for this model is 

£^ = Tr 1 1 VHS^'Cns)!?^^ - (^MV')H°7''(n5)V'} ; (42e) 
we do not include a mass term because, since 175 anti-commutes with ij^, the expression 

tph-i^i-l^il) (42f) 

is anti-self-adjoint, and vanishes inside Tr. It is easy to check that under the transformations 

P : x°) -> r]pij°tP{-x, x°) , 

C : ip —>■ r]ci>^'^ , 

T : 7P{x, x°) ^ VtA^P{x, -x°) , (42g) 

together with the gluon sector transformations discussed above, the Lagrangian density JL,^ is P odd, 

C odd, and T even. Although naively leading to conserved source currents for the gauge gluons, the 
_ 5 

model JL^ has chiral anomalies in the usual complex canonical quantization, which suggests that it may 
also be inconsistent in the more general total trace Lagrangian dynamics. This question requires further 
study. 

* When G ^ G', the models do not have a C or a CPT symmetry; this does not contradict the usual CPT 
theorem because we do not make the locality assumption that the fields commute at spacelike separations. In 
[1] it is shown that when the Lagrangians of this section are reinterpreted as complex field theory Lagrangians 
through a decomposition of the fields into symplectic components, then one finds a C and a CPT symmetry 
even for G ^ G', in agreement with the usual CPT theorem. 
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(5) Up to this point, our entire discussion has dealt with field theories in fiat space-time. Since the total 
trace Lagrangians of this section are all Lorentz invariant, they can be generalized to curved space- 
time by the standard prescription of replacing the Minkowski metric by a general metric g^,^, ordinary 
derivatives d\ by covariant derivatives Va which commute with g^^, etc. When this is done, the 
source term for the gravitational field equation will be a total trace energy-momentum tensor T'^'', 
defined by computing the variation of the total trace action S under an infinitesimal metric variation 
9niy g^u + ^9nv, according to 

S^S + 5S, ^^ = \ j d^x[g{x)Y'^T:^''{x)5g^,{x) . (43a) 

(Here g{x) is the negative of the determinant of the metric, and our notation follows Weinberg [24].) 
Standard arguments then show that T^" is covariantly conserved, 

VmT''" = > (43b) 

and, in the fiat space-time limit, the spatial integrals of the various components of T'*'' give the total 
trace Poincare generators, for example, 

H = y" d^xT°^ . (43c) 

We conclude that quaternionic field theories described by total trace actions can be consistently coupled 
to classical gravitation, but the total trace structure of the gravitational source term differs from that 
customarily assumed in the semi-classical theory of gravitation. 

(6) It is interesting to ask whether the Lagrangian L of Eqs. (30f) and (34b, c) has fermionic symme- 
tries, constructed in analogy with the fermionic symmetry of supersymmetric Yang-Mills theory [25,26]. 
This question can be investigated by studying the change of L under field variations parameterized by 
fermionic parameters, leading to a lengthy calculation, the results of which are compactly summarized 
as covariant divergence equations for the corresponding Noether currents. Let S^^-^ 2) be the fermionic 
currents 

<^(l,2) - 2 
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t ^ /t 

(1,2) -Q2' + ^(1,2) 



7°7^^ [7^,7.] 



(44a) 



which transform under operator gauge transformations as 



(44b) 



Let us define the covariant derivative Dx acting on a general operator O as 



D.O = d.O + B'O - OB, 



(44c) 



so that when O gauge transforms as 



O U'OU^ 



(44d) 



DxO transforms covariantly as 



DxO ^ U'{DxO)U^ 



(44e) 



Comparing Eq. (44c) with the definition of in Eq. (30c), we evidently have 



(DxO)^ = DxO^ . 



(44f) 



Forming the covariant divergence DxS^^^ 2) using Eq. (44f), we get the identity 



(1,2) - 2 



(I>aV'(1,2))^7^ + 7^ (^aV'Ji,2)) 



Q2 Q,2 



2) Q2 



+ 



Q,2 



^(1,2 



7°7^^ [7^,7.] • 



(45a) 



The first line of Eq. (45a) can be simplified by using the Dirac equations of Eq. (35), expressed in the 
form 

(£'aV'(i,2))^°7^ = V';i,2)™7° , (45b) 



to give 



4i 



(1,2) Q2 ^ QI2 "^(1,2) 



'^(1,2 



The second line of Eq. (45a) can be rearranged by substituting the identity 



(45c) 



(45d) 
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The contribution from s^^" vanishes by virtue of the Bianchi identities 



DxF^^ + D^Fx^, + D^F^x = , D'^F'^, + D'^F'y^^ + D'^F',^ = 



(45e) 



The contribution from the Kronecker delta terms is 



(1,2) Q2 



Q,2 -^(1,2) 



(45f) 



which can be simpHfied using the gauge field equations of Eq. (35) to give 



7%. 



(45g) 



Thus, putting everything together, we have 



1 / 4. F^^'^ F'f^^ 1 \ 1 
'Da^(i,2) = 2 y{i,2)-Q2 + -Qf^ ^(1,2) j "^^^2 f^"'^"^ 

+ (Wi,2)^^-^"'V';i,2))7V, 



(46a) 



and we see that even when the fermion mass m vanishes, the fermionic currents S^^ 2) s^'"® not covariantly 
conserved. Suppose, however, that there is either an operator gauge, or an asymptotic limit, in which 
the fermion fields V'(i,2) have the standard canonical anti-commutators of complex fields. In such a 
situation, we see from Eq. (35) that we would have J^^ = J'^, and the second line of Eq. (46a) would 
reduce to the singular commutator 

[V';i,2)'^17°7m, (46b) 

which vanishes in dimensional regularization. This argument suggests that Eq. (46a), despite the pres- 
ence of the gluon source current terms, may nonetheless have useful content. 

Let us now return to our main theme of total trace operator dynamics, and construct the total trace 
Hamiltonian form of the dynamics following from the scalar field Lagrangian of Eqs. (30e,f) and the fermion 
field Lagrangian of Eqs. (34b, c). From Eq. (30e), we get 



(5L 



sE 



--Q2 ^oe 



6L 

Pb' = — — 
' 6B' 



"qI2 ^oe ' 



(47a) 
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and so the total trace Hamiltonian density becomes 

3 



n = Tr 



Substituting 



• = Do(l)-Bo(l) + (t>B'o , 



G'2 



-C. 



(47b) 



(47c) 



Be = Foe + DiBq , B'^ — F^f + D[B'q , 
forming the total trace Hamiltonian H and doing a spatial integration by parts, we get 



with 



_ = j (fxTr 



3 



(48a) 



3 3 3 

E ^B,? - ^ E (^^™)' - ^0(^0 + E ^^PB,) 



£,m=l 



(48b) 



d^xTr 



? E (P'^'f - ^ E (^i-)' - B'oiX + E 



£=1 



<,m = l 



^=1 



and with Jo, Jq the components of the boson source currents Jv, Jl given in Eq. (32). Proceeding similarly 
in the fermion case, from Eqs. (34b, c) we get 



and so the total trace Hamiltonian density becomes 



(49a) 



1 3 ■ 1 ^ 

V'(2)V'(i) - V'(i)V'(2) - E ^oeBe - E ^o€-B; 



£=1 



Substituting 



V'(l,2) = £'oV'(l,2) - SoV'(l,2) + V'(l,2)-Bo , 



(49b) 



(49c) 



together with the second line of Eq. (47c), forming the total trace Hamiltonian and doing a spatial integration 
by parts, we now get 

H = H^(i 2) + Hb,Jo + ^B'-j^ , (50a) 
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with 

^=1 (50b) 
-V'(i)7V-D^V'(2) - (-D^V'(2))V7V(i) -m V'(2)nV(i) - ^(i)«7V(2) |- 

In Eq. (50a), Hb,Jo H^'.-j-^ are still given by Eq. (48b) with the substitution Jq —Jq, but now Jo, 

are the components of the fermion source currents Jv, Jl given in Eq. (35). 

The total trace Hamiltonian dynamics for '0(1)) and ■0(2) now takes the form of Eq. (5f), with no 

further complications. For the gauge potentials iJ^ and iJ^, however, we encounter the familiar problem 

that we are dealing with a constrained system, and so the canonical momenta are not independent. Focusing 

henceforth on the potential (the treatment of B'^ is completely analogous), we have a primary constraint 

Pbo = ^00 = , (51a) 

which is satisfied as an identity without use of the equations of motion. Differentiating Eq. (51a) with respect 
to time, we get the secondary constraint 

b ^ 
= Pbo = -j^ Hs,^„ = Jo + ^ biPB, , (51b) 

which is the same as the constraint arising from the Lagrangian equations of motion. Further time differen- 
tiation of Eq. (51b) leads to no further secondary constraints, since the equation 

^0 + ^ (j2 ^^^'s^) = (51c) 
can be rearranged [with use of Eq. (51b) and properties of the covariant derivative into the form 

D^J^. = , (51d) 

which is satisfied by virtue of the Lagrangian equations of motion for (j) or ^(1,2)- The constraint structure 
is thus completely analogous to that of a conventional Yang-Mills gauge field, for which the simplest way to 
realize a Hamiltonian dynamics is to use axial gauge [27,28], in which B3 is taken to vanish, 

B3 = . (52a) 
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In this gauge we have 



-^03 = ^o-Bs — d^Bo + [Bq, Bs] = —daBo , 
D3F03 = dsFoa + [B3, F03] = -dlBo , 

Fi3 = -^3^1 , F23 = -93^2 , (52b) 



and so the constraint 



Jo = J2^^ i-G^PBe) = ^i^oe = DiFoi + D2F02 - d^Bo (52c) 
e=i 1=1 

can be directly integrated to yield Bq and F03, giving (with 2:3 = z) 

Bo = -- / dz'\z - z'\{Jo + DiPB, + D2PB,)z' , 

^ J-00 

1 r°° z — zi 
Fo3 = i:G^ dz' {Jo + D^pb, + D2PB,).' ■ (52d) 

^ i-oo \Z-Z\ 

Substituting Eqs. (51b) and (52b) back into Hb^Jq, we get 



H 



= f d^xTr \ Yl (P^^)' - ^ (^03)^ - ^ {diB2 - d2B, + [BuB2]f 

^ [ £=1,2 



(52e) 



-^(%Bl)2-^(53S2)^^ 



with Fo3 given by Eq. (52d), and so only B\, B2 remain as independent dynamical degrees of freedom, with 
the corresponding independent canonical momenta Pbi,PB2- It is now completely straightforward to verify 
that the operator equations of motion obtained from the total trace Hamiltonian Hb^Jq of Eq. (52e), 

— — Co-Dl , —J — Co-t>2 , 

OPBi OPB2 

-sbT = "^^^ ' = "^^^ ' ^''^ 

are identical to the operator equations of motion obtained from the total trace Lagrangian. So we have 
achieved a consistent Hamiltonian dynamics. The generalized Poisson bracket of Eq. (6a), in axial gauge, 
now contains variational derivatives only with respect to the gluon variables Bi^2 and 3- We have not 
verified the Poincare generator algebra, but just as with the verification of the Hamiltonian equations of 
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motion which we have described, this should be a straightforward analog of the conventional Yang-Mills 
axial gauge calculation. 

At no point in the discussion have canonical commutation relations been used to get the operator 
equations of motion. They have been replaced in total trace dynamics by the constraints 

3 3 

= Jo + ^ DePBe , = ± Jo' + ^ DePB'^ , (54) 

and are a quaternionic field theory generalization, for independent left and right gaugings, of the constraints 
of Eq. (14b) of Sec. 2. We conjecture that any operator realization of Eq. (54) gives, via the total trace 
Hamiltonian formalism of Eqs. (47a)-(53), a consistent quaternionic field dynamics. 

5. Discussion 

In the preceding three sections, we have presented a generalization of standard quantum mechanics, 
providing a framework within which one can formulate quaternionic quantum field theories. We suggest in 
the final chapter of [1] that such theories may play a role in physics between the GUTS scale and the Planck 
mass, possibly providing a dynamics for preon models of quarks and leptons, and that the observed complex 
quantum field theories associated with the standard model and its grand unification are effective theories 
describing the asymptotic dynamics of the underlying quaternionic fields. 

These speculations aside, the concepts of total trace Lagrangian and Hamiltonian dynamics, and oper- 
ator gauge invariance, provide rich new possibilities for the formulation and study of quantum field systems. 
Among the many open questions which remain are: (1) Can the generalized bracket of Eq. (6a) be proved 
to satisfy a Jacobi identity? (2) What is the analog of the Dirac theory of constrained systems for operator 
gauge invariant systems? (3) What is the analog of BRST theory for such systems? (4) For the gauge models 
discussed in axial gauge in Sec. 4, what form does the total trace Hamiltonian take in non-canonical gauges, 
such as transverse gauge? (5) Does total trace dynamics, in the quaternionic case, correspond to a unitary 
time development, or is this only a property of complex quantum mechanics? (6) Can one usefully charac- 
terize the solutions of the operator constraints of Eq. (54)? (7) Given two scalar or fermion quaternionic 
operator fields, is there a criterion for determining whether they are related by a bi-unitary operator gauge 
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transformation? (8) Can one find a functional integration form of total trace dynamics, analogous to the 
Feynman path integral in complex quantum field theory? (9) More generally, what are useful calculational 
techniques for the theories formulated in this paper? 
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